Non-divergent pseudo-potential treatment of spin-polarized fermions under ID and 3D harmonic 
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Atom-atom scattering of bosonic one-dimensional (ID) atoms has been modeled successfully using a zero- 
range delta-function potential, while that of bosonic 3D atoms has been modeled successfully using Fermi- 
Huang's regularized .s-wave pseudo-potential. Here, we derive the eigenenergies of two spin-polarized ID 
fermions under external harmonic confinement interacting through a zero-range potential, which only acts on 
odd-parity wave functions, analytically. We also present a divergent-free zero-range potential treatment of two 
spin-polarized 3D fermions under harmonic confinement. Our pseudo-potential treatments are verified through 
numerical calculations for short-range model potentials. 

PACS numbers: 34.50.-s,34.10.+x 
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I. INTRODUCTION 

Recently, atom-atom scattering has received renewed in- 
terest since the properties of ultracold atomic (bosonic or 
fermionic) gases depend predominantly on a single atom- 
atom scattering parameter yj. This is the i-wave scattering 
length flj for a three-dimensional (3D) Bose gas |2l (or for a 
3D Fermi gas consisting of atoms with "spin-up" and "spin- 
down"), and the p-wave scattering volume Vp |3, 4] for a 3D 
spin-polarized Fermi gas. For a ID or quasi-lD gas, it is 
the ID scattering length a id fj, "ct], which characterizes the 
even-parity and odd-parity spatial wave function applicable to 
bosons and to spin-polarized fermions, respectively. In many 
instances, atom-atom scattering processes can be conveniently 
modeled through a shape-independent pseudo-potential |7, 8], 
whose coupling strength is chosen such that it reproduces the 
scattering properties of the full shape-dependent 3D or ID 
atom-atom potential. 

Fermi-Huang's regularized pseudo-potential (^ [T^ urn 
supports a single bound state for positive a^ and no bound 
state for negative Os- It has been used frequently to describe 
3D .s-wave scattering between two bosons or two fermions 
with different generalized spin. Busch et al. [12], e.g., de- 
rive the eigenenergies for two atoms under harmonic con- 
finement interacting through Fermi Huang's pseudo-potential 
analytically. Using an energy-dependent scattering length 
as{E), their results can be applied successfully to situations 
where Us is l arge and positive, i.e., near a Feshbach reso- 
nance I13lfl4l[l5ll . Building on these results, Borca et al. iTol 
use a simple two-atom model to explain many aspects of 
an experiment that produces molecules from a sea of cold 
atoms using magnetic field ramps |17]. In addition to these 
two-body applications, Fermi-Huang's 3D i-wave pseudo- 
potential plays a key role in developing (effective) many -body 
theories. 

This paper determines the eigenspectrum of two spin- 
polarized 3D fermions interacting through a regularized p- 
wave zero-range potential, parameterized through a single pa- 
rameter, i.e., the p-wave scattering volume Vp, under har- 
monic confinement analytically. Since wave functions with 
relative angular momentum / greater than zero have vanishing 
amplitude at r = (where r denotes the distance between the 



two atoms), our zero-range p-wave potential contains deriva- 
tive operators. Furthermore, it contains, following ideas sug- 
gested by Huang and Yang in 1957 1 11], a so-called regular- 
ization operator, which eliminates divergencies at r = that 
would arise otherwise. We show that our pseudo-potential 
imposes a boundary condition on the wave function at r = 
(see also Ref. ]18]); this boundary condition serves as an al- 
ternative representation of the p-wave pseudo-potential. Ear- 
lier studies, in contrast, impose a boundary condition at finite 
r, corresponding to a finite-range pseudo-potential with two 
parameters \ 19, 20]. The validity of our pseudo-potential is 
demonstrated by comparing the eigenenergies determined an- 
alytically for two particles under harmonic confinement with 
those determined numerically for shape-dependent atom-atom 
potentials. 



Due to significant advancements in trapping and cool- 
ing, to date cold atomic gases cannot only be trapped in 
3D geometries but also in quasi-2D and quasi- ID geome- 
tries I2ll I22L I23II . In the quasi- ID regime, the transverse mo- 
tion is "frozen out" so that the behaviors of atomic gases are 
dominated by the longitudinal motion. Quasi- ID gases can 
hence often be treated within a ID model, where the atoms 
are restricted to a line. To model ID atom-atom interac- 
tions, for which the spatial wave function has even parity, 
delta-function contact interactions have been used success- 
fully. In contrast to the 3D 5-wave delta-function potential, 
which requires a regularization, the ID delta-function pseudo- 
potential is non-divergent ]24]. To treat spin-polarized ID 
fermions, a pseudo-potential that acts on spatial wave func- 
tions with odd parity is needed. Here, we use such a pseudo- 
potential to determine the eigenenergies of two spin-polarized 
ID fermions under harmonic confinement analytically. Com- 
parison with eigenenergies determined numerically for shape- 
dependent ID atom-atom potentials illustrates the applicabil- 
ity of our ID pseudo- pote ntial. Our results confirm the Fermi- 
Bose duahty i25ll26l l27l l2i. l29ll in ID for two atoms under 
harmonic confinement. 



II. TWO INTERACTING ID PARTICLES UNDER 
HARMONIC CONFINEMENT 

Consider two ID atoms with mass m and coordinates zi and 
Z2, respectively, under external harmonic confinement, 



Vtrap{zi,Z2) 



-mcofizi+zi), 



(1) 



where (O^ denotes the angular frequency. After separating the 
center of mass and the relative motion, the Schrodinger equa- 
tion for the relative degree of freedom z, where z = zi—z\, 
reads 
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where 



ii\D 



l^idz^ 



-Viz) 



'- 2 2 



(2) 



(3) 



Here, V{z) denotes the ID atom-atom interaction potential, 
and /u the reduced mass, /j — m/2. 

Section In Al reviews the pseudo-potential treatment of two 
ID particles with even-parity eigenstates, i.e., two bosons 
or two fermions with opposite spin, under harmonic con- 
finement. Section lirBl determines the relative eigenenergies 
ZijQ for two spin-polarized ID fermions interacting through 
a momentum-dependent zero-range potential under harmonic 
confinement analytically. Section ITTcl benchmarks our treat- 
ment of the momentum-dependent zero-range potential by 
comparing with numerical results obtained for a short-range 
model potential. 



For \E^j^\ -^ oo (and correspondingly negative g^^), Eq. Q 
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reduces to lowest order to 



F+ - - 



ff 



2/j(fl 



+ \2' 

id) 



(6) 



which coincides with the exact binding energy of the pseudo- 
potential K^g„f/o(z) without confining potential. In Eq. (|6jl, 
fl[£, denotes the ID even-parity scattering length. 
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which is related to the ID coupling constant g^^ through 
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In Eq. (Q, k denotes the relative ID wave vector, k = 
^I/jEsc/K and E^c the ID scattering energy. The phase 
shift d^p is obtained by matching the free-space scattering 
solution for positive z to sin(fe + 5^^). The dashed line in 
Fig.QIa) shows the binding energy of the even-parity pseudo- 
potential without confinement, Eq. (|6|l, while the dash-dotted 
line shows the expansion of Eq. (|5Jl to next higher order. 

In addition to the ID eigenenergies E^j^, the eigen func- 
tions V|/]'q(z) can be determined analytically, resulting in the 
logarithmic derivative 



dz 






(9) 



A. Review of pseudo-potential treatment: Even parity 

The relative eigenenergies E^j^ corresponding to states 
with even parity (in the following referred to as even-parity 
eigenenergies) of two ID particles interacting through the 
zero-range pseudo-potential ^^^^^^^(z), where 



have been determined by Busch et al. iT 
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In Eq. (0), 5'''(z) denotes the usual ID delta function. The 
transcendental equation (|5} allows the coupling strength g^j^ 
for a given energy E^j^ to be determined readily. Vice versa, 
for a given ^j'^, the even-parity eigenenergies E^j^ can be 
determined semi-analytically. Figure ^a) shows the result- 
ing eigenenergies E^j^ of two ID bosons or two ID fermions 
with opposite spin as a function of the coupling strength g^. 
As expected, for vanishing interaction strength {g^^ — 0), 
the relative energies Ej^ coincide with the harmonic oscil- 
lator eigenenergies is,"^'' with even parity, E'"'^^ = (2n + ^)h(i>y, 
where n = 0, 1,- • •. 



This boundary condition is an alternative representation of the 
even-parity pseudo-potential V+ ^^^^^(z). 



B. Analytical pseudo-potential treatment: Odd parity 

Following the derivation of the even-parity eigenenergies 
by Busch et al. |12:], we now derive an analogous expres- 
sion for the odd-parity eigenenergies E^j^ using the zero-range 
pseudo-potential y^-.^„^^(z). 



y;seud.M) = ^^^zgw^^^'\^) 



dz 
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This pseudo-potential leads to discontinuous eigenfunctions 
with continuous derivatives at z = 0. We show that the loga- 
rithmic derivative of V(/[^(z) is well-behaved for z ^ 0+. In 
Eq. ( I lot , the first derivative acts to the left and the second to 
the right. 
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dz dz 



with the short-hand notation 
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Since V'^.^jj^^y (z) acts only on wave functions with odd par- 
ity (and not on those with even parity), we refer to V^^^^^jg{z) 
as odd-parity pseudo-potential; however, V_gjj^y(z) itself has 
even parity. Similar pseudo-potentials have recently also been 
used by others lElEJ Isoll. 

To start with, we expand the discontinuous odd-parity 
eigenfunction V(/]"o(z) in continuous ID odd-parity harmonic 
oscillator eigenfunctions (|)n(z). 
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where the c„ denote expansion coefficients, and 
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where a- = ^Jh/{^i(i>^). In Eq. (tlH . the L„ ' ' {z^ /al) denote 
associated Laguerre polynomials and the (|)„ (z) are normalized 
to one. 



\^„{z)\^dz^\. 
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The corresponding odd-parity harmonic oscillator eigenener- 

gies are 
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where « = 0, 1 , • • •. Inserting expansion ( I13t into Eq. (|2ji, mul- 
tiplying from the left with (j)*/(z), and integrating over z, re- 
sults in 
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The coefficients c„i are hence of the form 

dv:Ao) 
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dz 
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where the constant A is independent of n' . Inserting this ex- 
pression for the c„ into Eq. ( I17t leads to 
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If we define a non-integer quantum number v through 

and use expression \\A\ for the (|)„ (z), Eq. jl9> can be rewritten 
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where the z ^ 0+ limit is well-behaved. Equation (I21> can be 
evaluated using the identity 



-'"''>^>.r,-v)4-v,|.4 



(22) 



and the known small z behavior of the hypergeometric func- 

tionC/(-v,|,^)mj, 
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Using Eqs. J22l i and j23> in Eq. J21> . evaluating the derivative 
with respect to z, and then taking the z — > 0+ limit, results in 
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Replacing the non-integer quantum number v [see Eq. (I20H 
by E^jy/2h(a, — 3/4, we obtain the transcendental equation 
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which allows the ID odd-parity eigenenergies E^j^ to be de- 
termined for a given interaction strength g^j^. 

Solid lines in Fig.^b) show the ID odd-parity eigenener- 
gies /ij^, Eq. M5\ . as a function of ^j^. For g^j^ — 0, the 
eigenenergies E^j^ coincide with the odd-parity harmonic os- 
cillator eigenenergies E'"'^, Eq. ( I16> : they increase for positive 
gjQ ("repulsive interactions"), and decrease for negative gj"^ 
("attractive interactions"). 

Expansion of Eq. M5\ to lowest order for large and negative 
eigenenergy (implying positive g\jj), \E^jj\ -^ °°, results in 
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where the ID scattering length Oj^ is defined analogously to 
'Ad [with the superscript "+" in Eq. replaced by the su- 
perscript "— "]. The ID scattering length o^^ is related to the 
ID coupling strength ^j^ through 
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(27) 



The energy given by Eq. \26\ coincides with the binding en- 
ergy of the ID pseudo-potential K^gH^/f/z) without the confin- 
ing potential. A dashed line in Fig.fQb) shows iij^, Eq. M6\ . 
while a dash-dotted line shows the expansion of Eq. (I25> in- 
cluding the next order term. 

In addition to the eigenenergies Ej^, we calculate the eigen- 
functions 1|/]"q, 
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Figure 1 : Solid lines in panel (a) show the relative even-parity en- 
ergies £j'q [Eq. Js}] calculated using the pseudo-potential V'Jj^yjp(z) 

as a function of gj^. Solid lines in panel (b) show the relative odd- 
parity energies Ej^ [Eq. <25H calculated using the pseudo-potential 
^ttseudo^^) as a function of g^j^. Horizontal solid lines indicate the 
harmonic oscillator eigenenergies [with even parity in panel (a), and 
with odd parity in panel (b)]. Horizontal dotted lines indicate the 
asymptotic value of the eigenenergies E^j^ and E^j^ for gj'^ -^ ±oo 
and gj"^ — > ±<x), respectively. Dashed lines show the binding ener- 
gies E^j^, Eq. J5J, in panel (a) and fij^, Eq. J26> . in panel (b) of 
the pseudo-potentials V+^,,|^^(z) and ^pseudo<^' respectively, with- 
out confinement. Dash-dotted lines show the expansion of Eq. Js} 
[panel (a)] and Eq. <25> [panel (b)] including the next order term. 



Following steps similar to those outlined above, the logarith- 
mic derivative at z ^- 0+ reduces to 
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Figure 2: Relative odd-parity eigenenergies £[q for two particles un- 
der ID harmonic confinement as a function of the well depth d. Solid 
lines show the eigenenergies obtained by solving the ID Schrodinger 
equation, Eq. J2j, for the Hamiltonian given in Eq. J3} numerically 
using a short-range model potential, Eq. <31> . for a series of well 
depths d. Symbols show the eigenenergies obtained for the pseudo- 
potential ^^(,„^„(z)> taking the energy-dependence of the ID cou- 
pling constant g^^ into account, gj"^ = g'^j^{Esc) (see text). 



by Eq. Q with y[z) = V^gujoi^) [and vice versa, odd-parity 
energies fij^ can be obtained by solving the ID Schrodinger 
equation with V(z) — ^pseudiX^)^- ^^^ analytical treatment of 
two ID particles under external confinement thus confirms the 
Fermi-Bose duality for two ID particles under harmonic con- 
finement 121121123, mm. 



C. Comparison with shape-dependent ID atom-atom potential 

To benchmark the applicability of the odd-parity pseudo- 
potential V^^^^^^^^{z) to two ID atoms under harmonic con- 
finement, we solve the ID Schrodinger equation, Eq. (|2ji, for 
the Hamiltonian given by Eq. Q numerically for the shape- 
dependent Morse potential Vmorse{z), 



Vmorseiz) = de 
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(31) 



Our numerical calculations are performed for a fixed range 
parameter zq, zo = 11.65a.u., and for a = 0.35a.u.; these 
parameters roughly approximate the 3D Rb2 triplet poten- 



Equation (|55) is an alternative representation of the ID odd- tial |32]. The angular trapping frequency cOj. is fixed at 



parity pseudo-potential V 



pseudo V 



The even-parity eigenenergies E^j^ [Eq. (|5}] and the odd- 
parity eigenenergies E^j^ [Eq. J25H, as well as the logarith- 
mic derivatives [Eqs. Q and J29H are identical if the coupling 
constants of V+^^^^/z) and V^^.g„,/o(z) are chosen as follows. 
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This implies that even-parity energies E^j^ can be obtained by 
solving the ID Schrodinger equation, Eq. (|2ji, for H\d given 



10^'^a.u. (27tv, = CO,), and the atom mass m at that of the ^^Rb 
atom, implying an oscillator length a, of 1 12.5a.u., and hence 
a fairly tightly trapped atom pair. To investigate potentials 
with different ID scattering properties, we choose depth pa- 
rameters d for which the ID Morse potential supports between 
zero and two ID odd-parity bound states. Solid lines in Fig.|2] 
show the resulting ID odd-parity eigenenergies /Sj^ obtained 
numerically as a function of d. The corresponding eigenstates 
have "gas-like character", that is, these states would corre- 
spond to continuum states if the confining potential was ab- 
sent. 



To compare the odd-parity eigenenergies obtained numer- 
ically for the Morse potential V,„(,rse{z) with those obtained 
for the odd-parity pseudo-potential Vp^.g„^„(z), we follow 
Refs. [T?, 15]. We first perform scattering calculations for 
the ID Morse potential (no confinement) as a function of the 
relative scattering energy Esc for various depths d, which pro- 
vide, for a given d, the energy-dependent ID scattering length 
a^ij{Esc), where a^j^{Esc) = — tan(5["^(fe))/fe. Using the rela- 
tion between the ID scattering length a j^ and the ID coupling 
strength gj^, Eq. Ml\ . we then solve the transcendental equa- 
tion (I25> self-consistently for E^j^. 

Diamonds in Fig. |2 show the resulting odd-parity eigenen- 
ergies Eyjj for two ID particles under harmonic confinement 
interacting through the odd-parity energy-dependent pseudo- 
potential Vp,g„rf^(z) with g\jj = g\oiE>^c)- Excellent agree- 
ment between these eigenenergies and those obtained for the 
Morse potential (solid lines) is visible for all well depths d. 
We emphasize that this agreement depends crucially on the 
usage of energy-dependent ID coupling constants. In sum- 
mary. Fig. 12 illustrates that the odd-parity pseudo-potential 
y^seudo^-^) provides a good description of the eigenstates of 
two spin-polarized ID fermions under harmonic confinement 



for all interaction strengths, including g 
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III. TWO INTERACTING 3D PARTICLES UNDER 
HARMONIC CONFINEMENT 

Consider two 3D particles with mass m and coordinates n 
and r2, respectively, confined by the potential Vrrap(n ,^2), 
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(32) 



where ii>ho denotes the angular trapping frequency of the har- 
monic 3D confinement. The corresponding Schrodinger equa- 
tion decouples into a center of mass part, whose solution can 
be readily written down, and into a relative part. 
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Here, r denotes the relative coordinate vector {r — r2 — r\ ), 
V{f) the atom-atom interaction potential, and//||,' the 3D har- 
monic oscillator Hamiltonian, 
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The corresponding Schrodinger equation for the relative coor- 
dinate reads 



H3D^\>3D(f) = E^D^ioif) 



(35) 



Section lTlI Al brieflv reviews Fermi Huang's regularized s- 
wave pseudo-potential, while Section lTlI BI solves Eq. J35t for 
a regularized p-wave zero-range potential analytically. To il- 
lustrate the applicability of this p-wave pseudo-potential. Sec- 
tion IIII CI compares the resulting relative eigenenergies E30 
for two particles under harmonic confinement with those ob- 
tained numerically for a shape-dependent short-range model 
potential. 



A. Review of 3D pseudo-potential treatment: j-wave 

Using Fermi-Huang's regularized i-wave (l — 0) pseudo- 
potential y;=0^,(r)i,[IIl, 
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where 5'^'(r) denotes the radial component of the 3D 5- 
function. 
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and flj the 3D i-wave scattering length, Busch et al. \l'£\ de- 
rive a transcendental equation for the relative 3D eigenener- 
gies £30, 
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Here, fl/,„ denotes the oscillator length, fl/j„ — ■\/fi,/(/ja)/,„). 
Solid lines in Fig. |3ja) show the .s-wave energies E^o as a 
function of a^. For large and negative Et,d (and hence positive 
fly), an expansion of Eq. ( I38t to lowest order results in 
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which corresponds to the binding energy of V'^^^ jjjr) without 
the confining potential. A dashed line in Fig. [2 shows the 
energy given by Eq. (I39> . while a dash-dotted line shows the 
expansion of Eq. (13 8> including the next higher order term. 

Since only i-wave wave functions have a non-vanishing 
amplitude at r — Q, Fermi-Huang's regularized pseudo- 
potential leads exclusively to .s-wave scattering (no other par- 
tial waves are scattered). Equation ( I38> hence applies to two 
ultracold bosons under external confinement, for which higher 
even partial waves, such as d- or g-waves, are negligible. 

Recall that the irregular solution with / ~ diverges as 
r^'. The so-called regularization operator ^'" of the pseudo- 
potential y'Leudo^)^ ^l- * ^36l . cures this divergence. The solu- 
tions \|/3D (?) of two particles under external confinement obey 
the boundary condition 
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this boundary condition is an alternative representation of 
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B. Analytical 3D pseudo-potential treatment: p-wave 

The importance of angle-dependent p-wave interactions 
has recently been demonstrated experimentally for two potas- 
sium atoms in the vicinity of a magnetic field-dependent 



p-wave Feshbach resonance p31. Here, we use a p-wave 
pseudo-potential to model isotropic atom-atom interactions; 
treatment of anisotropic interactions is beyond the scope of 
this paper. 

We use the following p-wave pseudo-potential V'j^lj^^{r), 



<=L('^)=^rV,5(^)(r)Vri^r2, 



(41) 



where the coupling strength ^i "summarizes" the scattering 
properties of the original shape-dependent atom-atom interac- 
tion potential I34ll35ll. 
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Here, Vp denotes the /7-wave scattering volume 0], 

tan(5p(^)) 



Vn = lim - 



k^ 



(42) 



(43) 



hp the p-wave phase shift, and k the relative 3D collision 
wave vector Similarly to the ID odd-parity pseudo-potential 
ypseudcM)^ the first gradient V^ with respect to the relative vec- 
tor r acts to the left, while the second one acts to the right. 
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Just as the s-wave pseudo-potential V' ^^^^^^(r) does not cou- 
ple to partial waves with / ^ 0, the p-wave pseudo-potential 
^pseudoi'') '^"^^ ^^^ couple to partial waves with / 7^ 1 I3ql . 
Pseudo-potentials of the form ^i^Vy5^^-*(r)Vy^ have been 
used by a number of researchers before 13413^ BflBTT : dis- 
crepancies regarding the proper value of the coefficient gi, 
however, exist (see, e.g., Ref. |36]). Here, we introduce the 

regularization operator j^'"^ [Eq. ( I41H . which eliminates 
divergencies that would arise otherwise from the irregular p- 
wave solution (which diverges as r^^). A similar regular- 
ization operator has been proposed by Huang and Yang in 
1957 1 11 1; they, however, use it in conjunction with a cou- 
pling parameter gi different from that given by Eq. ( I42> . 
By comparing with numerical results for a shape-dependent 
model potential, we show that the pseudo-potential V''^J,j/„('^) 
describes the scattering behaviors of two spin-aligned 3D 
fermions properly (see Sec. lIIICV 

To determine the relative eigenenergies Et,d of two spin- 
polarized 3D fermions under harmonic confinement analyti- 
cally, we expand the 3D wave function x|/3£)(r) for fixed angu- 
lar momentum, / = 1, in continuous harmonic oscillator eigen 
functions (])„/„,, (r). 



mill 



(45) 



where the c„„,, denote expansion coefficients. The (])„/„,, (r) de- 
pend on the principal quantum number «, the angular momen- 
tum quantum number /, and the projection quantum number 



mi. 



and 



m^,^^,l„,{r)=E^^'Un,{^ 



pose 



2n + l+- ] hdiho, 



(46) 



(47) 



where n = 0, 1,- • •; / = 0, 1,- • • ,n— 1; and m/ = 0,±1,- • ■ ,±l. 
The isfninii i^) '^^^ ^^ written in spherical coordinates [r = 

('•,^,cp)], 



(48) 



^,dm,{r) = V47l/?„,(r) F,,„,(iJ,(p), 



where the F/m, ("djCp) denote spherical harmonics and the 
Rni{r) are given by 



Rnlir) = 




2' 



(2/ + l)!!V^Li'+'/2'(0)fli, 



r^ \ (/+1/2) /r2 



exp 



2«L 



'■ho 



(49) 



with 



(2/+l)!! = l-3-...-(2/+l). (50) 

The normalizations of R„i (r) and F/,,,, (i3-, (p) are chosen as 
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F;m,(i^,9)r sm^d'&d(? = l 



and 



\R„,{r)\\^dr^-. 
47t 



(51) 



(52) 



If we plug expansion ( 145 1 into the 3D Schrodinger equa- 
tion, Eq. ( I35> . for the Hamiltonian given by Eq. (I33> with 
Vi'') = K^eidoi^)^ multiply from the left with c^*,,,„;(r) [with 
/ = 1], and integrate over r, we obtain an expression for the 
coefficients c„/,„/ , 



Cn'm',{K'[ - E3d) = 



where 



v4^^k^L^,™;^«/M 



(53) 



Jr^O 



In deriving Eq. ( I53> . we use that 

V?[/?„,(r)Fz„,(iJ,(p)] = 
[V?-R„,(r)]F,,„,(i3,(p) +R„iir) [Vpy/„,,(i^,(p)] 



(54) 



(55) 



where the second term on the right-hand side goes to zero in 
the r — > limit. Since the gradients V^ in Eq. ( I53> act on 



arguments that depend solely on r, we can replace them by 
Cy-^, (where e,- denotes the unit vector in the r-direction), 



If we insert expansion (I63> into Eq. (I62> . evaluate the deriva- 
tives, and take the r ^ limit, we find 



(^n'm'i (K't ^ ^io) - 



-gl 



^KAO) 



1 3M 2 '^ 

2 3^3 I '" L^'^nm'Rntir) 



(56) 



Jr^O 



Equation ( I56> implies that the coefficients c„/,„/ are of the form 



C^,' ,J — A 



n'r ' 



(57) 



where A is a constant independent of n' . Plugging Eq. (I57> 
into Eq. ( I56> results in an implicit expression for the 3D ener- 
gies E^D, 









1 



(58) 



To simplify the infinite sum over n, we use expression ( I49> for 
the /?„/ (r), and introduce a non-integer quantum number v, 



h&ho a^ 



ho 



1 



r(-v) 



gl 



/Sr(-v-i)r(f)- 



(64) 



Using Eqs. ( I42> and ( I59> . we obtain our final expression for 
the relative eigenenergies Ejo for / = 1, 



V„ 



^ho 



n- 




'\) 


8r(- 


ElD 

2hm,„ 


+ 1) 



(65) 



Solid lines in Fig.|3lb) show the relative 3D eigenenergies 
Et,d, Eq. ( I65t . for two spin-polarized fermions under exter- 
nal harmonic confinement interacting through the zero-range 
pseudo-potential Vi^Xdo^) ^^ ^ function of the 3D scattering 
volume Vp. For vanishing coupling strength g\ (or equiva- 
lently, for Vp ~ 0), Ej,o coincides with the / = 1 harmonic 
oscillator eigenenergy. As Vp increases [decreases], £3/) in- 
creases [decreases]. 

Expansion of Eq. ( I65> for a large and negative eigenenergy 
(and hence negative Vp), \Et,d\ -^ °°, results in 



£30= 2v + / 



3V7t- 



2dr 



exp 



Using the identity 



r(3/2), 



For I = 1 , we obtain 



2a 



ho 



n=0 



-/a 



h&hf 



d'^'Hr^alS 



(59) 



n-V 



to, 



'ho aha. 



^ "ho 



gl 



Km 



(61) 



£■30^ 



(66) 



2A'(Vp)2/3' 

which agrees with the binding energy of V'^JjjJr) without the 
confinement potential. A dashed line in Fig.lSb) shows this 
binding energy, while a dash-dotted line shows the expansion 
of Eq. J65t including the next higher order. Compared to the 
eigenenergy of the system without confinement, Eq. ( I66t . the 
lowest eigenenergy given by Eq. ( I65> is downshifted. This 
downshift is somewhat counterintuitive, and contrary to the 
i-wave case. 

In addition to the eigenergies £30 of two atoms with / ~ 1 
under harmonic confinement, we determine the corresponding 
eigenfunctions \|/3D(r), 



the infinite sum in Eq. J60> can be rewritten. 



r(-v) 



1 33 



2aiu 



5 r' 
2 flL 



r^O 



gl 



m) 



where the r ^ limit is, as discussed above, due to the 
regularization operator of V^^^^^i^^{r) well behaved. Expres- 



sion \62l can be evaluated using the known small r behavior 
of the hypergeome trie function t/(— V, §, t-) 1311 . 



V 2 al 



Tt 



1 



r 

aho 



(2V + 3) 



r(-v) 



r(- 



)m 



+ 0{r). 



(63) 



V3Z)(?) ' 



r(-v) r 

[ahoT^^ "ho 



exp 



Ho 



u 



^ho 



:(67) 



which lead to the well-behaved boundary condition 



a3 



'1..2 



^Mr)) 



r^W3D{r) 



V„ 



(68) 



This boundary condition is an alternative representation of the 
pseudo-potential V'jJ^^^^{r), and depends on only one param- 
eter, that is, the scattering volume Vp. This is in contrast to 
earlier work 1 19, 20], which treated a boundary condition sim- 
ilar to Eq. ( I68> but evaluated the left hand side at a finite value 
of r, i.e., at r — re. The boundary condition containing the 
finite parameter r^ cannot be mapped to a zero-range pseudo- 
potential. References |38, 39, 40] discuss alternative deriva- 
tions and representations of boundary condition ( I68> . 
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Figure 3: Solid lines in panel (a) show the relative i-wave ener- 
gies £30 [Eq. <38H calculated using the pseudo-potential %j^^^^^(r) 
as a function of the scattering length Aj. Solid lines in panel (b) 
show the relative p-wave energies £3£) [Eq. <65H calculated using 
the pseudo-potential V'^,',^„ {j) as a function of the scattering volume 
Vp. Horizontal solid lines indicate the harmonic oscillator eigenen- 
ergies [for / = in panel (a), and for / = 1 in panel (b)]. Horizontal 
dotted lines indicate the asymptotic eigenenergies £3/5 [for as -^ ±0° 
in panel (a), and for Vp -^ ±00 in panel (b)]. Dashed lines show the 
binding energies, Eq. <39> in panel (a) and Eq. <66> in panel (b), of 
the pseudo-potentials VJ^'^^doi^) ™'l ^pTeudoi^)' respectively, with- 
out confinement. Dash-dotted lines show the expansion of Eq. <38> 
[panel (a)] and Eq. <651 [panel (b)] including the next order term. 



C. Comparison with shape-dependent 3D atom-atom potential 

To benchmark our p-wave pseudo-potential treatment of 
two spin-polarized 3D fermions under harmonic confinement, 
we solve the 3D Schrodinger equation, Eq. J35> . for the 
Hamiltonian given by Eq. ( I33> numerically for the shape- 
dependent Morse potential Vi„urse{f), Eq. (13 U with z replaced 
by r and zo replaced by ro. As in Sec. Ill CI our numerical 
calculations are performed for ro = 11. 65a. u., a — 0.35a.u., 
K)/!o = lO^^a.u. (2nVho — W/io), and m ~ m(^^Rb). The well 
depth d is chosen such that the 3D Morse potential supports 
between zero and two / = 1 bound states. Solid lines in Fig.|3 
show the resulting 3D eigenenergies £30 with / = 1 obtained 
numerically as a function of the depth d. 

To compare the / = 1 eigenenergies obtained numerically 
for the Morse potential Vmorsei'') with those obtained for the 
p-wave pseudo-potential V''^J,^/„('^), we follow the procedure 




500 1000 1500 2000 
potential depth d/hi^ho 

Figure 4: Relative 3D eigenenergies £"30 with 1=1 for two spin- 
polarized fermions under 3D harmonic confinement as a function 
of the well depth d. Solid lines show the eigenenergies obtained 
by solving the 3D Schrodinger equation, Eq. <35> , for the Hamilto- 
nian given in Eq. <33> numerically for a short-range model poten- 
tial, Eq. )31t with z replaced by r and zo replaced by ro, for a series 
of well depths d. Symbols show the eigenenergies obtained for the 
pseudo-potential V j'2tdg{r), taking the energy-dependence of the 3D 



scattering volume Vp into account, Vp 



-VpiEs, 



(see text). 



outlined in Sec. Ill CI that is, we first determine the energy- 



dependent free-space scattering volume Vp{Esc), Vp{Esc) = 
— tsin{5 p{k))/k^, for the 3D Morse potential (no confinement) 
as a function of the relative scattering energy E^c for various 
well depths d. We then solve the transcendental equation ( I65l l 
self-consistently for E^d. Diamonds in Fig. show the re- 
sulting / = 1 eigenenergies Ej,o for two 3D particles under 
harmonic confinement interacting through the / = 1 energy- 
dependent pseudo-potential V'^J,^,/'^) with Vp = Vp (£■«■). Ex- 
cellent agreement between these eigenenergies and those ob- 
tained for the Morse potential (solid lines) is visible for all 
well depths d. We emphasize that this agreement depends 
crucially on the usage of energy-dependent 3D scattering vol- 
umes. Figure 13 illustrates that the p-wave pseudo-potential 
^psJudoi^) describes p-wave scattering processes properly. 



IV. SUMMARY 

We determined the eigenspectrum for two ID particles un- 
der harmonic confinement interacting through a momentum- 
dependent zero-range potential. This pseudo-potential acts 
only on states with odd-parity, and is hence applicable to 
the scattering between two spin-polarized ID fermions. We 
showed that a basis set expansion in continuous functions 
can be used to determine the eigenenergies and discontinu- 
ous eigenfunctions of two ID particles under harmonic con- 
finement interacting through the odd-parity pseudo-potential 
^meudoi^)- ^^^ divergence-free treatment confirms the Fermi- 
Bose duality in ID for two particles. 

We also determined an implicit expression for the eigenen- 
ergies £30, Eq. J65l l. and eigenfunctions \\i3D{r), Eq. J67l i. 
of two spin-polarized 3D fermions under harmonic confine- 
ment interacting through a momentum-dependent zero-range 



9 



potential. Similar to studies of two atoms with / = fl3'. 14t 
Il5. .Idl . our analytical expressions might be useful in under- 
standing the behavior of two confined spin-aligned fermions, 
including physics near Feshbach resonances. The p-wave 
pseudo-potential used in our study contains derivative oper- 
ators as well as a regularization operator; the former is needed 
to construct a true zero-range potential (since 1 = 1 solutions 
go to zero as r approaches zero, see above) while the latter is 
needed to eliminate divergencies of the irregular /?-wave solu- 
tion (which diverges as r^^). We showed that our zero-range 
potential V'^j',jo('^) imposes a boundary condition at r = 0, 
Eq. (I68> . which depends on a single atomic physics parame- 
ter, that is, the scattering volume Vp. This boundary condition 
is an alternative representation of Vt:J^^^^^{r). 

Similarly to Fermi-Huang's regularized i-wave pseudo- 
potential, the /7-wave pseudo-potential used here might find 



applications in developing effective many-body theories for 
ultracold spin-polarized Fermi gases. Such theories will have 
to carefully investigate how to implement renormalization 
procedures needed in numerical calculations. 

Note added: After submission of this paper we became 
aware of a related study by Stock et al, see quant-ph/0405 153] 
which derives Eq. (I65> starting with a pseudo-potential ex- 
pressed as the limit of a 5-shell. 
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